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Abstract 

We use QCD sum rules to estimate the universal form factors describing the 

semileptonic B decays into excited charmed resonances, such as the 1~ and 

3- 

2~ states and D2 belonging to the ~ 2 ^^^^^ quark doublet, and the 

5- 

2~ and 3^ states D2' and belonging to the sf = - doublet. 
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1 Introduction 



In the heavy quark {Q = c, b) infinite mass hmit (mg —* oo) Quantum Chromodynamics 
exhibits symmetries that are not present in the finite mass theory: heavy quark spin and 
flavour symmetries as well as the velocity superselection rule 0. These approximate 
symmetries allow to organize the spectrum of physical states comprising one light anti- 
quark and one heavy quark in multiplets of definite parity P and total angular momentum 
Si of the light degrees of freedom. 

The lowest lying multiplet consists in the meson doublet with = | , corresponding 
to the vector 1~ and the pseudoscalar 0~ state. The doublet can be described by a 4 x 4 
Dirac matrix 

H = ^-^[p;!" - Pi,] (1) 

where v is the heavy meson velocity, P*^ and Pa are annihilation operators of the 1~ and 
0~ Qqa mesons (a = 1, 2, 3 for u, d and s); for charm, they are D* and D, respectively. 

The nearest mass multiplets are the sf = doublet, comprising the positive parity 
l"*" and 0"^ states, and the sf = I"*" doublet which includes the positive parity 1"*" and 2+ 
states. In the charm sector three of such states have been identified: the state 1^2(2460) 
is the narrow 2"*" meson with = l"*^; moreover, there are two l"*" mesons with masses 
m^o = (2422.2 ± 1.8) MeV § and m^*o = (2461^^^ ± 10 ± 32) MeV |]; they can 
be identified with members of the multiplets predicted by the Heavy Quark Effective 
Theory including some mixing between them. Evidence for such states has also been 
collected in the beauty sector |p. From the theoretical viewpoint these states have been 
the subject of intense scrutiny: the role of the \^ doublet (O"*", 1+) in some applications 
of chiral perturbation theory has been considered in |^ and in [0; their properties have 
been studied both by QCD sum rules |], |10] and quark models |jlT . 



In this letter we investigate some properties of the next heavy meson multiplets, the 
sf = I doublet including two mesons with = 1^ and 2", and the sf = | doublet 
which comprises the states with = 2~ and 3~. We estimate the universal form factors 
describing, in the infinite heavy quark mass limit, the semileptonic B decays into such 
multiplets, and consider the contribution of these processes to the inclusive semileptonic 
B decay width f\ 

^The operators in have dimension § since they contain a factor ^Jmp in their definition. 

review on the problems related to inclusive and exclusive semileptonic B decays can be found in 
ref.0. 



We follow the QCD sum rule approach |T^, which has been applied to similar problems 
in the past ||^, [1^, [T^ However, as discussed in the following, in the application of the 
method to high-spin states several difficulties appear in identifying the range of parameters 
needed in the sum rule analyses, due to the peculiar features of the considered states and 
of their interpolating currents. In order to overcome such difficulties, we make use of 
information coming from other theoretical approaches, namely constituent quark models 
predicting the heavy meson spectrum. The final result, although affected by a sizeable 
theoretical uncertainty, nevertheless is useful for assessing the role of high-spin meson 
doublets in constituting part of the charm inclusive semileptonic B decay width. 



2 Effective meson operators and quark currents 



The effective operators describing the sf = | and sf = | meson doublets are given 
respectively by 0: 
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(2) 
(3) 



where D* represent annihilation operators of the mesons with appropriate quantum num- 
bers. In order to implement the QCD sum rule programme, we need quark currents with 



non- vanishing projection on these states. They have been investigated in ref.|ll0[ and are 
given by the following expressions: 
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;§)-; J"" = 3- 



jaf3\ rpaf3X, fj.ucrj^ 
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(4) 
(5) 
g(6) 
(7) 



where D'^ is the covariant derivative: = — i g A'^ , and represents the 

transverse component of the four-vector with respect to the heavy quark velocity v. 



*For a review see ifTa 



Gi = — {G-v)v^ . The tensors T""^'^^ and T°'^^'^^'' are needed to symmetrize indices 
and are given by 



- \ [gfg^'^gr + gfg^i'' + i^^gT) , (9) 



with (7°^ = (7°^ — v°'v^ 



As discussed in in the mg oo hmit the currents in eqs.(^)-(0) have non- 

vanishing projection only to the corresponding states of the HQET, without mixing with 
states of the same quantum number but different si content. Therefore, we can define a 
set of one-particle-current couplings as follows: 
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= /2 e*-^ 
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= /2 e*"^ 


(12) 


- 
— 


(1)- 


JP 


= 3- : 


<d;{v, e) 


jaf3X 
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= /3 e*-^^ , 


(13) 



where e are the meson polarization tensors. The couplings fi are low-energy parameters, 
determined by the dynamics of the light degrees of freedom. Since the two pairs (/i, /2) 
and (/2, /a) are related by the spin symmetry, in the sequel we only consider /i and f2- 

3 Two-point function sum rules 

To evaluate the parameters /i and /2 let us consider the two-point correlators 

t J (fx e"^""-" < |T(jt"(x)J^(0)| > = ni{uj)gf (14) 

I /dSe-*-^-^<o|T(jt-'3(^)j>-(o)|o> = ii2{uj){g79t + grg^t' -\9f9r) 

(15) 

given in terms of IIi and 112, scalar functions of the variable uo. 

As extensively discussed in the literature, the QCD sum rule method amounts to eval- 
uate the correlators in two equivalent ways. On one side the Operator Product Expansion 



(OPE) is applied for negative values of u; the expansion produces an asymptotic series, 
whose leading term is the perturbative contribution (computed in HQET), followed by 
subleading terms parameterized by non perturbative quantities, such as the quark conden- 
sate: < qq >, the gluon condensate: < asG^yC^^ >, the mixed quark-gluon condensate, 
etc. On the other side, one evaluates the correlators by writing down dispersion relations 
(DR) for the scalar functions Ili{uj) and 112(0;); they get contributions by the hadronic 
states, in particular by the low-lying resonances with appropriate quantum numbers. To 
get rid of radial excitations and multiparticle states, one performs a Borel transform on 
both sides of the sum rule, which enhances the low mass contribution of the spectrum; 
moreover, assuming quark-hadron duality, one identifies, from some effective continuum 
threshold Uc, the hadronic side of the sum rules with the perturbative result obtained by 
the OPE. In the final sum rule, only the contributions from the physical to the continuum 
threshold appear: the low mass resonance on one side, the OPE truncated at Uc on the 
other. 

Applying the method to the correlators ([14D and ([l5| ) we get two borelized sum rules 
for the parameters /i and /2: 

/Je-/'' = A r-daa'e-'^ + (16) 

TT"^ JO lO 

f2^-A./E' ^ }^ da a'e-'^/^' . (17) 
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Here the parameters Ai and A2 are defined by the formulae: Ai = m/j* — and A2 = 
m£)*i —rric, rric being the charm quark mass; therefore, the parameters Ai and A2 represent 
the binding energy of the states Dl and D2', which is finite in the infinite heavy quark 
mass limit. On the other hand, ujic and u!2c represent the effective thresholds separating 
the low-lying resonances from the continuum; E and E' are parameters introduced by 
the Borel procedure. Relations for the mass parameters Ai and A2 can be obtained by 
taking derivatives of the sum rules (|T6D and ([T7|): 

1 fi^lc 



/ da a^e-"!^ 
Jo 



^Tdaa^e-^/'^ + < ^""^^ > 



Jo ' ' ' 32 
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rda aV^/-' 

A2 = ^u^^ . (19) 

daa'e--/^' 

Jo 

There is an important point deserving a discussion, and it concerns the high dimen- 
sionality of the interpolating currents J" and J°^, which has two consequences on the 



structure of the sum rules (]T^)-(]T7|) and ([TB|)-(|T^. First, the spectral functions in eqs.(|TBp- 
(|T^) and (jT^)-(|T^) have large powers, and therefore the perturbative contributions in the 
sum rules are very sensitive to the continuum thresholds ujic and 002 c- The second ef- 
fect consists in the absence of the contributions from low-dimensional condensates, which 
implies (neglecting high-dimensional condensates) complete duality between the perturba- 
tive and the hadronic contributions to the sum rules. Such two effects cannot be avoided 
in our analysis, and are typical of the sum rule approach to high spin states |jl6[. In 
our case they have the main consequence of not allowing to determine simultaneously the 
couplings fi and the mass parameters Aj, due to the critical dependence on the continuum 
thresholds. Therefore, we adopt the strategy of getting the values of the mass parameters 
from other determinations, and then to fix the thresholds from eqs.(p!8D-(p!9D and comput- 
ing fi from (p!6D-(p!7]). Admittedly, this is a hybrid procedure, which nevertheless allows 
us to estimate both the current-particle matrix elements and the universal semileptonic 
form factors, as discussed in the next Section. 

While experimental information on the = | and = | doublets are not available 



so far, there are studies concerning such states based on constituent quark models [17 



They suggest that the mass of the 3~ (cm) state D3 is rriDs = 2.83 GeV or = 2.76 
GeV, whereas the mass of the corresponding (bu) state is =6.11 GeV. Assuming a 
spin splitting of ~ 40 MeV in the charm sector, as suggested by the same models, we can 
give to the mass of the | state the value of 2.78 GeV, e.g. nearly 0.8 GeV above the 0~ 
doublet (the same value comes from the analysis of the beauty meson spectrum). This 
implies for the parameter A2 a value in the range A2 — [1.2 — 1.4] GeV, considering the 
determination of the analogous binding energy of B and D mesons [|T5|. As for Ai, we fix 
it to Ai ^ [1.2 — 1.4] GeV, according to similar considerations. 

Let us consider Ai and A2 related to the thresholds Ui and to the Borel parameters 
Ei by eqs.([T8|)-(|l9D. There is a range of Borel parameters and thresholds where the 
chosen binding energies can be obtained. In particular, while the dependence of Aj on 
the Borel parameters is quite mild, so that the range Ei = [1 — 1.5] GeV can be chosen, 
the dependence on the thresholds, as expected, is critical: one has to choose Ui in a quite 
narrow range [1.6 — 1.8] GeV to obtain Aj. However, this choice is not unappropriate, 
since it suggests that the continuum threshold is above the mass of the corresponding 
resonance by nearly the mass of one pion. 



After having fixed Aj and the ranges of Ei and of cUj, from eqs.(p!6|)- ([T7|) we can obtain 

5 ~ 7 

the values of the couplings ff. fi = [0.6 — 0.8] GeVa and /2 = [1.2 — 1.6] GeVa. Notice 
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that, at odds, e.g., with the leptonic constants related to the matrix elements of the 
quark axial currents on the 0~ state, the couplings fi do not have an immediate physical 
meaning, as they represent the projections of the interpolating currents on the orbitally 
excited meson states. Nevertheless, they play an important role in the determination of 
the form factors, as we discuss in the next Section. 



4 Universal form factors from three-point sum rules 



There are two universal form factors describing the semileptonic B decays into the excited 
negative pa 
the decays 



negative parity charmed resonances with = | and = | . The first one, ri, governs 



B DH 

B d;i 



(20) 
(21) 



in the heavy quark limit. The second one, T2, describes in the same limit the decays 

B D;'iiye (22) 
B ^ Dsiue . (23) 

It is straightforward to write down the semileptonic matrix elements for the transitions 
(|20|)-(p3D, by applying, e.g., the trace formalism |T^. One obtains: 



<Dl{v',e)\{V-Ar\B{v)> = ^niBniD* n{y) ^- [{e*v) [v^ - v"' + 



(24) 



(25) 



for the decays (^) and (0), while for the decays (p^) and (|23|) the relevant matrix 
elements can be written as: 



<D*^{v\e)\{V -AY\B{v)> 



-^niBmD*! T2{y) el^v' 



[2(1 -y') 



5 h ^ 



<D^{v\e)\{V-AY\B{v)> = ./W^^^ r^^y) e^^vl' g^^\l + y) - v\' ^ 



(26) 



(27) 



In these equations the weak current is {V — A)^ = C7'^(l — 75)&, y = v -v' and Ti{y), T2{y) 
are the universal form factors. 

At the zero- recoil point v = v' the matrix elements in (^)-(^) vanish, as expected 
by the heavy quark symmetry. As a matter of fact, for B decays into spin 2 and spin 
3 states, at least one index of the final meson polarization tensor is contracted by the 
B four-velocity v, and therefore the product vanishes for v = v'. The spin symmetry 
requirement being verified in the matrix elements, the Isgur-Wise form factors ri and T2 
are not required to vanish a.t v ■ v' = 1. 

One can attempt an estimate of the form factors ri 2 by three-point function sum rules, 
considering the correlators (relevant for the matrix elements (|2^) and (^)): 

J d^x d^z ^-i(^v-x-Jv'-z) ^ Q |T(jt«(^)(y _ AY{^)Jr,{x)\ > = 

= ie^'^^^Vfiv'^^^{uj, cj') + Si(cj, u')w'^v^ + ... (28) 
J d^x d^z ^-K^v-x-Jv'-z) ^ Q \T{J^^\z){y - AY{^)J^[x)\ > = 

= ie^^^Pvrv'piw'^gl + M;'^^°)fi2(^, ^') + w'^w^v^'^^i}^ . + •••(29) 

where = f " — J5 = qij^b; the dots represent other Lorentz structures which are 
not relevant for the subsequent analysis, since we only consider Qi and 0,2- 

Since the scalar functions Qj depend on two variables, one has to perform double DRs 
and double Borel transforms, which introduces, for each sum rule, two Borel parameters 
E and E'. The resulting equations read: 

n{y) = ^ J^^ e^/^^^^/^' r rdada'e-^'^-'^l^'hM a')e{a, a') (30) 
2\/27r2fiF •^'o Jo 



tM = -^47£^e^/^+^^/^' r rdada' e--/^--/^'/.2(a, a')e{a, a') (31) 
V27r2/2F -^0 Jo 



where 
hi{(T, a') 



1 



(|/2 - 1)3/2 



+ - 2yaa' a'{a + a') 



2{y-l) 3 



h2{a, a') = ^y^-^^^y2_ 1)5/2 [5^' - 3^V^(% - 1) + (2y2 - 2y + l){a'' + 3aa"] 



(32) 



and 



X(a, a') = 0(a2 + a'2-2yaa') , 
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(33) 



with 0(x) the step function. 

In eqs.(0) and (|3T|) the parameter A represents the mass difference between the low 
lying multiplet si = and the heavy quark. The integration region can be expressed 
in terms of the variables 



0"_ 



cr + cr' 



a -a' 



2 

and one can choose the triangular region defined by the bounds: 



0< cr+ <coiy) (34) 




V y+ 1 V 2/ + 1 

As to the upper limit in the integration interval for uj^ we adopt 

for the two cases studied in this letter (we use, according to the two-point sum rule 
analysis Ud = ujc2)- 

We use the value F = 0.21 GeV^/^ which is obtained by QCD sum rules ^ with 
as = (the same order which we consider in the present analysis). Moreover, we use 
A = 0.5 GeV, with the threshold in the B channel Uc = 0.7 GeV. As for the charm 
channel, we use Uic = oj^c = 1-6 — 1.8 GeV. 

We can now numerically determine the form factors r^, using the above equations. 
The result for the universal function Ti{y), obtained within the uncertainties discussed 
above, is that this function, in the whole kinematical region relevant for the decays (pOl)- 
(|2T|), is less than 10~^, which implies that, in the infinite heavy quark mass limit, the 
semileptonic B transitions into the S£ = | doublet have a very small decay width. The 
situation is different for the universal function T2{y), which is depicted in fig.^ where the 
shaded region corresponds to the results obtained by varying the parameters A, A2, ujc 
and uj2c in the ranges quoted above. The form factor T2, at the zero recoil point y = 1, 
is in the range T2(1) = 0.10 — 0.20, with a mild dependence that can be neglected, 
within the accuracy of the sum-rule method. Although it is difficult to reliably assess the 
theoretical accuracy of this result, it is interesting to observe that a form factor in the 
range quoted above implies that the semileptonic channel is experimentally accessible. 
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Figure 1: Universal form factor T2{y) 



5 Semileptonic decay rates 



Using the parameterization of the B matrix elements in eqs.(p6D and ( pTl) we can work 
out the expressions of the widths of the decay modes ( ^2]) and (^31), which are respectively 
given by: 



dT 
dy 



dV 



{B^ D 



— {B^Dsii^e) 
dy 



7207r3 



7207r3 



[r2(y))^(l/-l)^(l/+l)5[(l+r^)(7i/-3)-2r(%2_3y+3)] 

(37) 

{r2{y)f{y - 1)^ (y + 1)^ [(1 + r2)(ll + 32/) - 2r{lly + 3)] 

(38) 



with r = — ^. Using m^g = 2.78 GeV, mo" = 2.74 GeV and T2{y) = 0.15, we get 



rriB 



T{B D*'iui) ^ T{B D^iue) ~ 4 x IQ-^^ GeV 



(39) 



and 



B(B D 



2 '^'^e) 



B(B D 



1 X 10" 



(40) 
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Therefore, although small, semileptonic B decays to the | doublet are within the reach 
of the running B factories, and could be experimentally observed, since the final mesons, 
as discussed in the next Section, are expected to be rather narrow. 

As for B decays to the | doublet, due to the small value of the universal function ri, 
the semileptonic widths turn out to be negligible at the leading order in the ^ expansion 
(a discussion of the role of next-to-leading corrections for semileptonic B decays to excited 
charm mesons can be found in [0). 



6 Remarks on strong decays of orbitally excited charm 
states 

One might expect that the states in the multiplets sf = | and | , being significantly 
higher in mass than the low- lying sf = | multiplet, are rather broad. However this 
should be only true for the sf = | states. As a matter of fact, the = 2" and 
= 1~ states can decay into the 0~ or 1~ heavy meson plus one pion by P— wave 



3 



transitions, which implies a kinematical suppression of the order of , where A^ ^ 1 

GeV is the typical chiral symmetry breaking scale. Taking into account that, for the 
charmed mesons, — 0.68 GeV, we expect a kinematical phase space suppression, for 
this decay channel, of ~ 0.3. 

On the other hand, for the mesons belonging to the multiplet = | , the decay into 
the low lying heavy meson and one pion occurs by F— wave transitions: the kinematical 

\Ptt\'^ 

suppression is ~ , which numerically means a reducing factor ^ 0.07. Since the decay 

mode with one pion in the final state is expected to dominate the decay width, one may 
guess that the 3~ and 2~ mesons belonging to the = | doublet are rather narrow Q. 

To render these conclusions more quantitative, let us consider the effective lagrangian 
describing, in the chiral effective theory for heavy mesons the strong couplings of the 
multiplet H'^'^ to the pion and the multiplet H: 

C = ^Tr {hH^" [h{D^, D,}Ax + k2 {D^D^A, + D^D^A^)] 7^5} + h.c. (41) 

X 



^Similar conclusions are reached in \l{ 
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where H is the = | multiplet containing the and 1 low-lying states, and ki^2 
effective couplings; moreover 

Ax = \[ed^i~id>^e] (42) 
and ^ = exp[zy^]. Putting k = ki + k2, one obtains for the two-body decay widths: 

r(D,^D.) = (43) 

TiD,^D*n) = ^^M^a^ (44) 

T{D;'^D7v) = (45) 

r ^2 ^ D n) = 46 

57r/^A4 niD*: 

with ^ 132 MeV. The value of k is unknown; however, on the basis of QCD sum rule 
results for similar couplings one may assume k G [0.25, 0.5]. In correspondence to the 
lower bound in this range we get 

r(D3 ^ {D, D*)tt) ~ 32 MeV (47) 
T{D;' D*n) ~ 15 MeV (48) 

where we have assumed the mass splitting of 40 MeV between the 3^ and the 2^ mesons 
in the multiplet. There are other decay channels contributing to the full widths, but 
the corresponding partial widths are expected to be much smaller: for the decay modes 
with one pion and an excited positive parity D resonance in the final state, occuring by 
D— wave transitions, we estimate a width of 1-2 MeV; for the decay modes with two pions 
and a heavy meson in the final state we expect, in the infinite heavy quark mass limit, a 
negligible contribution. 

We can therefore conclude that reasonable estimates for the full widths of the S£ = | 
resonances are as follows: 

r(D3) = 35 - 140 MeV (49) 
T{D;') = 17 -70 MeV , (50) 



a consideration which suggests the presence of a not too broad peak in the Dir and D*tt 
channel in the region of 2.8 GeV. 
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This conclusion, together with the result of a branching fraction of semileptonic B 
decays to the | doublet of the order of 10~^, encourages the experimental investigation 
at the currently running i?-factories as well as at the hadronic facilities. 
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